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Parton energy loss in Au+Au collisions at RHIC energies is studied by numerically solving the 
relativistic Boltzmann equation for the partons including 2^2 and 2^2 + final state radiation 
collision processes. Final particle spectra are obtained using two hadronization models; the Lund 
string fragmentation and independent fragmentation models. Recent, preliminary tv° transverse 
momentum distributions from central Au+Au collisions at RHIC are reproduced using gluon-gluon 
scattering cross sections of 5 — 12 mb, depending upon the hadronization model. Comparisons with 
the HIJING jet quenching algorithm are made. 
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Exciting, preliminary ir° transverse momentum distri- 
butions from Au+Au collisions at ^snn — 130 GeV at 
RHIC have been recently reported When normal- 

ized to the mean number of binary collisions, the ratio 
of the central to the peripheral distributions increases up 
to px ~ 2 GeV/c, reaches a maximum value that is less 
than 1 and then decreases for p? > 2 GeV. This prelim- 
inary data reveals a suppression of the production of 7r°s 
in central Au+Au collisions for px > 2 GeV/c. 

The preliminary 7r° transverse momentum distribution 
has been reproduced by several calculations that included 
jet quenching (3|Q|. Jet quenching occurs when a high 
energy jet passes through a medium and looses energy 
from the induced non-abelian radiation. Recently, in- 
tense theoretical activity has been devoted to calculating 
the energy loss of a fast parton traveling through QCD 
media |^-|^]. It has been shown that transverse momen- 
tum distributions of hadrons at large px are sensitive to 
the total energy loss of the fast partons || . Jet quenching 
has been proposed as one of the signals of the formation 
of a quark-gluon plasma (QGP) P-full. Recent pQCD 
motivated calculations that include the nuclear effects of 
shadowing and energy loss via a modified fragmentation 
function have shown that the tt° data can be reproduced 
with a constant energy loss of dE/dz = 0.25 GeV/fm [||. 
Using a similar approach, but with a dE/dz that depends 
upon the number of collisions, it was shown that the data 
can also be reproduced with n — 3 — 4 average number 
of scatterings 

While there have been many theoretical studies of jet 
quenching, few simulation models have incorporated this 
physics. The HIJING simulation model (l2) and the two 
other models which use the HIJING algorithms, such as 
HIJING/BB ||] and AMPT @ incorporate a simple jet 
quenc hing mechanism. The HIJING jet quenching algo- 
rithm [p.2| assumes a simple gluon splitting scheme with a 



fixed energy loss dE/dz. The energy loss for gluon jets is 
twice that of the quark jets. A jet can only interact with 
locally comoving matter (strings) in the transverse di- 
rection and the points of the interactions are determined 
by the probability dP(l) — dl/\e~ l / x , where A is the 
mean free path and I is the distance the jet has traveled 
between collisions. When an interaction occurs (with a 
string), the medium induced radiation is simulated by 
forming a gluon kink in the string with AE — IdE/dz of 
the jets energy. A jet can interact with the surrounding 
medium until it exits the system or its pt fall below a 
certain pq cut-off. As the energy of a nuclear collision in- 
creases, more partons are produced, and this jet quench- 
ing algorithm leads to the production of many additional 
low energy gluons (gluon kinks). In a default HIJING 
jet quenching calculation with dE/dz = 2 GeV/fm, the 
charged hadron multiplicity per participant was shown to 
increase with the colliding energy much faster than the 
data (l5Hl|]. 

The effects of the parton scattering phase where only 
elastic parton interactions are included have been studied 
with several models; ZPC Q, AMPT @(l|L MPC |(J 
and GROMIT Q. Calculations revealed |9| that elas- 
tic parton interactions only slightly decrease the Et and 
multiplicity and produce very little elliptic flow when typ- 
ical gluon-gluon pQCD cross section of a gg ~ 3 mb and 
dN g /dy\ v= o ~ 200 are used. A large elliptic flow is only 
obtained when gross cross sections or large gluon densi- 
ties are used p0|] . Inelastic processes, like gg <-> ggg, have 
been shown to be be important |22]-|2~4|| for thermalizing 
the partons and in the radiative energy loss of the high 
energy jets. 

In this work, a model is introduced that describes 
the parton scattering phase for high energy heavy ion 
collisions by numerically solving the Boltzmann equa- 
tion for a distribution of partons with 2^2 and 
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2^2 + final state radiation interactions. The 2 — + 
2 + final state radiation processes are needed to effec- 
tively simulate jet quenching and the dynamical scat- 
terings of the partons lead to a reduction of the Et at 
mid-rapidity. As a result, the incident energy dependence 
of the multiplicity per participant in this approach is sim- 
ilar to HIJING without jet quenching and is consistent 
with data 

The initial parton distributions in our model are ob- 
tained from the HIJING event generator which produces 
around 190 gluons at mid-rapidity for a typical Au+Au 
collision at ^Jsmn — 130 GeV. The nuclear shadowing 
effects in HIJING are included in all calculations in this 
paper. The system of partons is then evolved in time 
using the relativistic Boltzmann equation, 



p»dJ(x,p)=C, 



(1) 



where f{x,p) is the distribution function of the partons 
and C is the collision integral. For the collision integral, 
2^2 and 2^2 + final state radiation processes are 
included. The distribution function, f(x,p) is assumed 
to be the sum of the particles, 
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where N tes t is the number of "test particles" and A is the 
actual number of partons produced in a collision. The 
momentum of the partons is determined by HIJING and 
the space-time coordinates are calculated using simple 
uncertainty relations. The formation time for partons is 
taken to be a Lorentzian distribution with a half width 
t = E/mlj,, where E and tut are the parton energy and 
transverse mass, respectively [[19). 

Monte-Carlo simulation is used to solve the Boltzmann 
equation. In this simulation, the on-shell partons are 
evolved in time along straight lines between collisions, 
where all of the collisions are time-ordered in a global 
frame. After each collision, the flavor, position and mo- 
mentum of the outgoing partons are determined and the 
list of possible collisions is updated. The space-time evo- 
lution of the system continues until there are no possible 
collisions between the partons. 

The "parallel-ensemble" method has been widely used 
to simulate the low energy nucleus- nucleus collisions p5| . 
In this method, collisions are determined using the "clos- 
est distance approach", where a collision occurs if the 
minimum relative distance b re i for any pair of particles 
in their center of mass frame becomes less than inter- 
action range as given by yj a/n. Here, a is the total 
parton-parton cross section. However, this method vio- 
lates Lorentz invariance due to the nature of action at 
a distance; a problem that has been studied by sev- 
eral authors (l^j2(],|2l],^6| . Two solutions to this prob- 
lem have been proposed. The "full-ensemble" method 
replaces each particle by Nt es t particles that interact 
with a reduced cross section a /Ntest p0|,p7 28 1. In the 



limit of Nt es t — * oo, one obtains the locality in config- 
uration space. Another solution is the "local-ensemble" 
method where the probability for one pair of the 

test particles to collide during the time interval of At in 
the small volume element AV is given by 



W = avAt/(N t est^V). 



(3) 



Here v is the relative velocity of the scattering particles. 
In the limit of AV — > 0, At -> 0, N test -> oo, the so- 
lutions will converge to the exact solutions of the Boltz- 
mann equation. In this study, the full-ensemble method 
is used with AT test = 6. 

The parton-parton interactions included in this simu- 
lation are 



q + q' ^q + q', g + g <-> g + g, 
g + g^q + q, g + q^ g + q. 



(4) 



The cross sections for these processes are given by lead- 
ing order (LO) perturbative QCD (LOpQCD) |j and 
explicitly take into account the running coupling con- 
stant a s (Q 2 ). Here, the Q 2 scale is chosen to be the 
squared transverse momentum transfer of the scattering 
process Q 2 = p\ and a s (Q 2 ) is evaluated at the scale of 
1 GeV 2 , when pr < 1 GeV. Since these cross sections 
diverge when pt — * 0, a Debye screening mass mjj = 0.6 
GeV |M and quark thermal mass m„ = 0.2 GeV are in- 
troduced to regulate the propagators |2(]]24|] . While these 
values should change with time, they are taken to be 
constant for simplicity. In this paper, two sets of results 
are obtained by multiplying these cross sections with two 
different factors (K- factors); K — 1.0 and K = 2.5. For 
example, a if-factor of 1.0 yields a gg ~ 5.3 mb for the 
gluon-gluon cross section, a gq ~ 2.0 mb for the quark- 
gluon, and a qq ~ 0.5 mb for the quark-quark, and a 
JT-factor of 2.5 yields 2.5 times the values for K = 1. 

The 2-^2 + final state radiation processes are mod- 
eled using PYTHIA algorithms |32j, where the two out- 
going partons are evolved with time-like branching tak- 
ing into account angular ordering. This approach has 
been used to study the gluon jet production in e + e~ 
interactions. During the branching, the life times of 
newly branched partons are obtained using the uncer- 
tainly principle, tf orrn ~ 1/Q, where Q represents the 
off-shellness of the parton in its rest frame |3^]. The 
minimum virtuality for the final state radiation is cho- 
sen to be 0.5 GeV |33|]. The number of emitted gluons 
(amount of final state radiation) is sensitive to this cut- 
off parameter. Using a Q 2 — 1 GeV 2 with K = 2.5 
yielded less final state radiation and hence less energy 
loss. A systematic study will be given later. The maxi- 
mum virtuality in parton-parton scattering is assumed to 
be Q 2 max = 4pj, where pt is the transverse momentum 
transfer of the scattering process. Off-shell partons are 
not allowed to interact with other partons in the system. 
While the final state (time-like) radiation is occurring, 
possible interactions with other partons are ignored, in 
contrast to the parton cascade model VNI E2] . 
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FIG. 1. Rapidity distributions of mini-jet gluons for 
Au + Au collisions (b < 4.48 fm) at ^sWn = 130 GeV. The 
solid line denotes the distribution from the initial condition 
(HIJING no jet quenching), and the dashed and dotted line 
represent the distributions after rescattering with a K factor 
of 1.0 and 2.5 respectively. The dotted-dashed line is the re- 
sult from HIJING default calculation {dE/dz = 2.0 GeV/fm). 

In Fig. [I], the rapidity distributions of the mini-jet glu- 
ons for central Au+Au collisions at ^/snn = 130 GeV 
are shown. The solid curve represents the initial rapidity 
distribution of gluons obtained from the HIJING model 
without jet quenching; dN/dy g _ y=0 = 190. With this 
initial distribution, the dynamical evolution of the par- 
tons with 2 — > n interactions enhances the gluon mul- 
tiplicity, where dN/dy gyV= o = 210 for K = 1.0, and 
dN/dy g ,y = o = 250 for K = 2.5. In comparison, HIJING 
with jet quenching with dE/dz = 2.0 GeV/fm yields 280 
gluons near mid-rapidity. As observed, these values are 
sensitive to the magnitude of the parton-parton cross sec- 
tions. 
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FIG. 2. Transverse energy rapidity distributions of par- 
tons for Au+Au collisions (b < 4.48 fm) at ^/snn = 130 GeV. 
The solid line denotes the distribution from the initial con- 
dition and the dashed^ = 1.0) and dotted (K — 2.5) line 
represents the distributions after rescattering. HIJING with 
default jet quenching result is shown by the dotted-dashed 
line. 
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FIG. 3. Rapidity correlation for parton-parton scattering 
for Au+Au collisions (b < 4.48 fm) at -/sJTn = 130 GeV. 
K — 2.5 is used. 

In Fig. ||, the rapidity distributions of mini-jet trans- 
verse energy dEr/dy in central Au+Au collisions at 
■J snn = 130 GeV is shown. After the partons are dy- 
namically evolved, the dEx/dy of the partons is reduced 
by approximately 13 GeV for K = 1.0 and 27 GeV for 
K = 2.5. The loss in Et results from the strong ra- 
pidity correlation in the parton-parton collisions. The 
rapidity correlation for the case of K — 2.5 is plotted 
in Fig. |^. As the partons with similar rapidities col- 
lide, the momentum is redistributed from transverse to 
longitudinal and the Et decreases. The degree of the 
correlation depends on the parton-parton cross sections, 
where the correlation becomes stronger as the cross sec- 
tions decrease. Most collisions in the simulation occur 
around the pr ~ 1-0 GeV/c. In comparison, in Fig. |^, 
the HIJING jet quenching scheme increases the dEx/dy. 
In the HIJING jet quenching scheme, the medium in- 
duced radiation is simulated by forming a collinear gluon 
(kink in the string) . In this approach, little of the trans- 
verse momenta is redistributed. 

Fig. |i| shows the gluon transverse momentum distri- 
butions from our calculations at the Snn = 130 GeV 
for central Au+Au collisions together with HIJING jet 
quenching calculations. The inclusion of 2 — * 2 + 
final state radiation processes results in the reduction 
of high pt partons and a slight increase in the low 
pr < 1 GeV region. This effect is sensitive to the size of 



3 



the parton-parton cross sections and the assumed energy 
loss dE/dz. 
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FIG. 4. Transverse momentum distributions of gluons 
for Au + Au collisions (b < 4.48 fm) at ,/sjvjv = 130 GeV. 
The solid line denotes the distribution from the initial condi- 
tion and the dashed and dotted lines represent the distribu- 
tions after rescattering with K = 1.0 and K = 2.5, respec- 
tively. HIJING results are also shown in squares (dE/dz = 2.0 
GeV/fm) and circles (dE/dz = 0.25 GeV/fm). 

Once the partons have finished interacting, two differ- 
ent hadronization models are used to obtain the hadron 
spectra; the Lund string fragmentation and independent 
fragmentation models. The results of these two models 
are compared with one another and with the 7r° data. 

When using the Lund string fragmentation model |32) , 
the string configurations are maintained throughout the 
evolution of the partonic system by determining the 
string configuration from the color amplitude of each 
parton-parton collision. In the PYTHIA manual |3l), de- 
tailed explanations of the possible string configurations 
for the relevant parton-parton scattering processes are 
provided. 

As shown in Fig. [| the HIJING without jet quench- 
ing calculation (solid line) overestimates the data at 
Pt > 2 GeV/c as consistent with other pQCD moti- 
vated parton model calculations (|||. J et quenching is 
not seen in the low pt region, 1 < pr < 2 GeV/c, as 
expected. As a result of the dynamical evolution of the 
partons, the data can be reproduced with the K = 2.5 
cross section set. In addition, HIJING with jet quench- 
ing dE/dz = 0.25 — 2.0 GeV/fm is able to reproduce 
the data. The dominant contribution of the ir° with 
Pt ~ 2—4 GeV/c comes from the hadronization of quarks 
with pt > 4 GeV/c. Thus, the ir° distribution is sensitive 
to the energy loss of the quarks and it is very important 
to include the appropriate qq and qg interactions. 
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FIG. 5. A comparison of neutral pion transverse momen- 
tum distributions for central Au + Au collisions (b < 4.48 fm) 
at ys^" = 130 GeV with the PHENIX data §| is presented. 
The solid line denotes the distribution from the initial con- 
dition and the dashed line represents the gluon distribution 
after rescattering with the K factor of 2.5. The HIJING re- 
sults for dE/dz = 2.0 GeV/fm (dotted line) and dE/dz = 0.25 
GeV/fm (dotted-dashed line) are also shown. 

Since explicit local color connections are not believed 
to be maintained in a high density partonic system, the 
independent fragmentation model is also used to com- 
pute the final hadron spectra. The independent frag- 
mentation model fragments all of the partons indepen- 
dently. As compared to the Lund string model, this 
model does not incorporate long range correlations be- 
tween the partons. Within this hadronization scheme, 
the high px 7r° data can be reproduced with K = 1.0 
as shown in Fig. |^. Not shown is a calculation with 
K = 2.5 which had a larger energy loss, underestimat- 
ing the data. Since the long-range correlations between 
partons are absent, the final hadron distribution strongly 
reflects the parton distribution. Changing the fragmen- 
tation model in HIJING from the default Lund string 
model to the independent fragmentation model, the HI- 
JING calculations with jet quenching show similar be- 
havior. Above pr > 2.0 GeV/c, calculations assuming 
dE/dz = 0.25 GeV/fm are able to reproduce the data. 
This value is the same as that obtained from the pQCD 
motivated parton model calculations 

The total hadron yields from our model near mid- 
rapidity are smaller than the HIJING without jet quench- 
ing values by approximately 4% when using Lund string 
fragmentation model and are larger by approximately 
2% when using the independent fragmentation model. 
We have checked that the model yields similar results at 
•J snn = 200 GeV. Our calculations are therefore con- 
sistent with data on the energy dependence of the 
hadron yield near mid-rapidity. 
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FIG. 6. A comparison of neutral pion transverse momen- 
tum distributions for central Au + Au collisions (b < 4.48 
fm) at v^snn = 130 GeV with the PHENIX data || is pre- 
sented. The solid line denotes the distribution from the ini- 
tial condition and the dashed line represents the distribution 
after rescattering using K = 1. The HIJING calculations 
(independent fragmentation model used) with jet quenching 
for dE/dz = 2.0 GeV/fm (dotted line) and dE/dz = 0.25 
GeV/fm (dotted-dashed line) are also shown. 

In this paper, we study parton energy loss in cen- 
tral Au+Au collisions at RHIC by numerically solv- 
ing the Boltzmann equation with 2 <-> 2 and 2 — > 
2 + final state radiation parton processes. Preliminary 
data on the ir° pr distributions at RHIC were fit with 
gluon-gluon cross sections of 5-12 mb, depending upon 
the hadronization model. This approach provides the 
maximum incoherent bound for energy loss in compari- 
son to the LPM effect. 

With the given HIJING initial condition, the parton 
scatterings reduce the parton cIEt / dy\ y= o by 5-10%. The 
magnitude of the change in the Et is sensitive to the scat- 
tering cross section. Interactions between particles that 
are strongly correlated in rapidity lead to a redistribution 
of the momenta from the transverse to the longitudinal 
direction and are important in lowering the Et and the fi- 
nal hadron multiplicity. Similar effects have been pointed 
out in regard to the suppressed production of open charm 
in the pre-equilibrium stage |34| . 

To test the sensitivity of the 7r° spectra to the 
non-perturbative hadronization process, two different 
hadronization schemes were used. Larger parton-parton 
scattering cross sections were needed to reproduce the 7r° 
distributions when using Lund string fragmentation than 
when using independent fragmentation. The spectra pro- 
duced by indepndent fragmentation are more sensitive 
to changes in the parton distributions, due to the lack 
of the long-range correlations introduced by the strings. 
These two different hadronization schemes therefore re- 
quire different amounts of parton energy loss in order to 
reproduce the data. 

While we have investigated the energy loss in the par- 



ton scattering phase, future calculations should also ad- 
dress the change in the pt distributions of the particles 
from the late hadron gas stage. The effect of the par- 
ton scattering phase using other initial conditions |35|-|39|l 
should also be explored. 
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